Stellar Populations 


For many modern applications, one is not concerned with the 
evolution of a single star, but with an entire set of stars. There 
are a number of sophisticated computer codes that track this, but 
the calculations are actually fairly straightforward (and largely 
analytic). 


In order to make the mathematics a bit more straightforward, 
let’s choose a fiducial mass, mı, which can be the mass of the 
Sun. Such a star will have a mean main-sequence luminosity, 
lı, and a main-sequence lifetime of Tı. We will then define the 
variable, m as the dimensionless mass of a star, i.e., for a star of 
mass M, m = M/m. 








First we need is an Initial Mass Function (IMF) of stars, which 
gives the number of stars born as a function of mass. The original 
IMF was that by Salpeter (1955) which was a simple power law, 


o(m)dm = ġım 7t?) dm = pım~ dm (31.01) 


with x = 1.35 (or s = 2.35), and limits between 0.1 Mọ and 
100 Mo. Other commonly used IMFs include that by Miller & 
Scalo (1979), Kroupa et al. (1993), and Chabrier (2003). (Most 
of the differences are at the low-mass end of the function.) The 
constant ġı is simply there for normalization purposes, so that a 
star cluster (or galaxy) with total mass Mo satifies 
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The second piece of information we need is the mass-luminosity 
relation for main-sequence stars. As we have seen, this is a 
roughly a power-law, with a slope of a ~ 3.88 at the faint end, 
and a slightly flatter relation at higher masses. For simplicity, 
we'll use a single power law connecting mass to luminosity on the 
main sequence 


la = lim (31.03) 


where the subscript d represents the luminosities of dwarfs (i.e., 
main sequence stars), and a ~ 3.5. 


The third item needed is the length of time stars spend on the 
main sequence. The lifetime of a star is simply proportional to 
the energy available to the star divided by rate at which the star 
is emitting its energy, i.e., its luminosity. Since the available 
energy is proportional to the available mass, and the luminosity 


is proportional to mass (though the mass-luminosity relation), 
m m 

(=) x m! = r=rnm!? (31.04) 
where 7, is the lifetime of the fiducial star with m = mı. Al- 


ternatively, this equation can be inverted. After a time t, the 
turnoff mass of a single-age stellar population will be 
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The final pieces of information that are needed are the lifetime 
of a typical post-main sequence (giant) star (Tg), the average 
luminosity of a giant star, £g, and the mass of a remnant star 
(typically a white dwarf, w). All of these can be derived (approx- 
imately) from models of stellar evolution. Note that Tz < Ta, so 
the stellar main-sequence lifetime is also roughly the star’s total 
lifetime. 





We will now consider the evolution of a set of stars all born (with 
the same metallicity) at the same time. This is called a Simple 
Stellar Population (SSP). The properties of more complicated 
systems can be inferred from the summation of several SSPs. 


Luminosity Evolution 


Let’s first calculate the luminosity evolution of a cluster (or 
galaxy) of stars with total mass Mo all born at the same time. 
The total luminosity of main sequence stars is easy to compute: 
all we have to do is sum up the luminosity of all stars on the 
main sequence. The lower end of the sum is the minimum mass 
of an energy-generating star; the upper end is defined by the 
main sequence turnoff. In other words, 
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Note that since the exponent is negative, and low mass stars 
live (essentially) forever, the last term in the above equation is 
negligible. Thus, the total luminosity of dwarf stars is 
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Calculating the total luminosity of giant stars is equally simple. 
First, note that the timescale for giant branch evolution is much 
faster than that for main sequence evolution. Thus, the key is 
to estimate the number of stars currently turning into giants; 
when this number is multiplied by the length of time a typical 
star remains a giant, and the mean luminosity of the star, the 
result is the total giant star luminosity. Now consider: the rate 
at which main sequence stars turn into giants is defined by how 
many stars are at the main sequence turnoff, and how much of the 
main sequence is eaten away per unit time interval. According 
to our definition of the initial mass function, the number of stars 
at the main-sequence turnoff is 
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and the number of stars turning into red giants during a time At 
is 
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so the total luminosity of giant stars is 
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If we now substitute for t for m using (31.05) take the derivative, 
and combine terms, the result is 
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Thus, the total luminosity of the stellar system, as a function of 
time, is 


br = Lp + bq = “Oe E) lil) E 
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(31.12) 


This can be simplified a bit if we define G(t) as the ratio of giant 
star luminosity to dwarf star luminosity, t.e., 


a(y = So = a aNlate (41) (31.13) 
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Since a > 1, the exponent of time in (31.13) is significantly less 
than one. Hence writing the expression for G(t) in this way 
illustrates that the variable is only a weak function of time. This 
expression can then be further simplied by using (31.03) and 
(31.05) 
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Note that the first part of this equation is very close to one. 
Furthermore, the term in brackets is simply the ratio of the total 
energy generated by the star on the giant branch to the total 
energy the star generated on the main sequence. Since stars turn 
off the main sequence when ~ 10% of their total fuel is gone, 
but eventually burn ~ 70% of their fuel, G(t) ~ 6. (From an 
observers point of view, it’s a bit more complicated, since most 





of the light a giant star produces is red, while the main sequence 
light may be blue. Thus, the exact value of G(t) depends on 
the wavelength of observation: in the blue, G(t) ~ 1, but in the 
red, G(t) > 10.) Using this notation, the luminosity of a stellar 
population, as a function of time, is 


Lr(t) = Lol) {1+ G(t)} 
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Figure 5. The relative 
contribution of the 
various evolutionary 
stages to the integrated 
bolometric luminosity of 
a simple stellar popula- 
tion as a function of age 
(lower scale) and turnoff 
mass (upper scale). The 
various evolutionary 
stages are defined in the 
text, composition and 
mass loss parameter are 
indicated, and s = 2.35. 




















Figure 4. The luminosity evolution of ` 
SSP's for three choices of the IMF | 
slope s, normalized at AM,,,=0fort - 
= 15 Gyrs. The lower line shows the — 
luminosity evolution of a population - 
with constant star formation rate, s = 
2.35, and upper IMF cutoff at 50 Mo. 


Photometric Mass-to-Light Ratio 


In addition to a population’s total luminosity, there are several 
other quantities which can be computed from first principles. 
First is the photometric mass-to-light ratio. Recall the total lu- 
minosity of main sequence stars in a stellar population is 





Lp Modif Ime =e (31.06) 
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where Min is the turnoff mass of the main sequence. As stated 

previously, a — x is usually greater than zero, so the last term of 

this equation is negligible; stars at the low end of the mass se- 

quence do not contribute significantly to the system’s luminosity. 


On the other hand, low mass stars can be an important part of 
the cluster’s mass. From (31.02), the visible mass of stars in a 
galaxy is given by 
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(neglecting the mass in invisible stellar remnants). If x > 1 
(as is normally assumed), the exponents in (31.17) are less than 
zero, and the value of M is dominated by the last term in the 
equation. Physically, this means that most of the stellar mass of 
a population resides in low mass stars. 


The combination of (31.06) and (31.17) means that it is virtually 
impossible to determine a population’s photometric mass-to-light 
ratio from observations alone. One can always drive up M/L by 
postulating the existence of low-mass stars which add mass (via 
31.17), without adding luminosity (31.06). 





Stellar Evolutionary Flux 


An interesting property associated with stellar populations is the 
stellar evolutionary flux, i.e., the number of stars passing through 
any (post main-sequence) phase of evolution at a given time. Just 
as the amount of water flowing down a river is controlled by the 
spillway of a dam, the number of stars evolving through various 
stages of evolution is controlled by the rate at which these stars 
move off the main sequence. (All other rates are much faster 
than this.) As we have seen previously, this rate is 
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From (31.05) 
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so the stellar evolutionary flux is 
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Naturally, the number of stars evolving through any phase of evo- 
lution in a galaxy is proportional to the number of stars in the 
galaxy. The best way to take this dependence out is to normal- 
ize the evolutionary flux to the size of the galaxy. This cannot 
be done with mass, since we do not know the total number of 
stars that are present in the system. However, we do know the 
population’s luminosity, which is defined in (31.16). So, when 
we normalize (31.20) to this luminosity, we arrive at the popula- 
tion’s luminosity-specific stellar evolutionary flux. After a 
bit of math, this quantity comes out to be 
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Note that the exponent of time is less than 1; thus B is rather 
insensitive to the age of the stellar population. (The difference 
between the stellar evolutionary flux of a 7 Gyr stellar population 
and a 12 Gyr population is less than a 25%.) Similarly, the 
specific stellar evolutionary flux is not very sensitive to the exact 
value of x. 
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Figure 2. The evolutionary flux b(t) 
from Eq. (3) for three choices of the 
IMF slope s, normalized at b = 1 fort = 
15 Gyrs. For each value of s the upper 
line gives the death rate, b(t - tp ys) 
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Detailed calculations show that the luminosity-specific stellar 
evolutionary flux is B ~ 2 x 107+} stars yr7! Le for all old (or 
moderately old) stellar populations. This makes is easy to pre- 
dict the number of any post-main sequence type of star present 
in a galaxy: since you know B, and you can measure the lu- 
minosity of the galaxy, all you need is the lifetime of the phase 
in question. For example, the lifetime of a planetary nebula is 
T ~ 25,000 years. If a galaxy is observed to have an absolute 
luminosity of £ = 10!! Lo, then the number of PN in the galaxy 
is N(PN) ~ B - £ - T = 50,000. 
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The plot above shows how a population’s bolometric-luminosity 
specific stellar evolutionary flux changes with the systems age 
and IMF. Note that the x-axis (population age) is a log quantity, 
while the y-axis is linear. All old (> 3 Gyr) stellar population 
evolve at rate of ~ 2 x 107! stars per year per solar luminosity, 
independent of age, IMF, or metallicity. 


[Renzini & Buzzoni 1986, in Spectral Evolution of Galaxies, p. 195] 


Two notes. First, the theoretical calculation of the luminosity- 
specific stellar evolutionary flux are based on bolometric luminos- 
ity. One rarely has access to the entire electromagnetic spectrum 
of a system of stars — most often the amount of luminosity sur- 
veyed in a star system is given in a specific filter, such as V or 
B. Fortunately, the conversion between V-luminosity and bolo- 
metric luminosity is very insensitive to the details of the stellar 
population. So, a bolometric correction of ~ —0.85 can convert 
V to bolometric magnitude. 
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Second: a common idea these days is that all (or many) plan- 
etary nebulae are formed out of binary star interactions. If so, 
then the previous calculation may not work, since it does not 
take into account binary-star interactions. (Interestingly, how- 
ever, the estimate does appear to come close to predicting PN 
numbers.) Similarly, if not all stars evolve through a given stage, 
then the calculation may be in error. 


Mass Loss from Stars 


Another useful quantity to know is the rate of mass loss from 
stars as a function of time. As before, the key to calculating 
this is to realize that almost all the mass lost from stars comes 
during the post main-sequence phase of evolution. Thus, the rate 
at which a population of stars loses mass is just by the rate at 
which stars move off the main sequence times the amount of mass 
each star loses. If m is the initial mass of the star, and w the 
mass of the stellar remnant (i.e., white dwarf), then the mass 
ejection rate is 








E(t) = Nin (Main — w) = Mob (Main) "(ren — w) (31.22) 


After performing the same substitutions as was done for the cal- 
culation of stellar evolutionary flux, this simplifies to 
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Once again, if we normalize to luminosity via (31.16), we can 
derive the luminosity specific mass loss rate 
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Numerically, this works out to a mass loss rate of ~ 0.02Mo 
per Gigayear per unit solar luminosity for an old (101° yr) stellar 
population. When integrated over the lifetime of a galaxy, the 
result is that ~ 15% of the initial stellar mass will be lost during 
a Hubble time. 
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More Complex Systems 


Real galaxies usually contain multiple stellar populations, each 
with its own mass, age, and metallicity. Computing the behavior 
of such a system simply means summing variations SSP compo- 
nents. It would, however, be useful to some zeroth order approx- 
imation to know the mix of stellar populations within a galaxy. 





A common way to do this is to simply assume that the star- 
formation rate of a galaxy declines exponentially with time, with 
some decay rate T, t.e., 


M(t) « e7*/7 (31.25) 


Irregular and late-type spiral galaxies are assigned large values 
of T, suggesting an almost constant star-formation rate over the 
history of the universe. Elliptical galaxies have small values of 
T, thus reflecting the fact that their star-formation has e-folded 
away long ago. Note however, that, many/most galaxies likely 
undergo a succession of discrete starbursts, with low-level star 
formation occuring in between times of extreme activity. So these 
models may be of only moderate use. 
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Problems in Photometric Evolution Calculations 


The above analytic solutions are only useful for guidance. To 
compare with observations, predictions must be made in specific 
bandpasses, or for specific absorption lines. This requires numer- 
ical calculations which include 


sets of stellar isochrones, detailing the precise number of stars at 
any position in the HR diagram as a function of age. 


sets of model stellar atmospheres, giving the emitted spectrum 
(or broadband color) of each star, as a function of its tempera- 
ture, surface gravity, and metallicity. 


an understanding of various non-stellar processes that may effect 
the emergent spectrum from galaxies. Such effects include inter- 
nal extinction due to dust, the presence of emission-lines from 
ionized gas, and (for high-redshift objects) absorption due to the 
intervening intergalactic medium (i.e., the Lya forest). 





Today, the most controversial parts of these codes include their 
treatment of thermally pulsing AGB stars, and the proper mod- 
elling of the horizontal branch and post-AGB stars. The former is 
important since it involves extremely luminous stars, which, due 
to their mass loss, may be obscured by circumstellar extinction. 
The latter is important for modeling old stellar systems which 
have very few blue stars. For these objects, a small change in 
the morphology of the horizontal branch (i.e., the ratio of red 
to blue objects), or the number of slow-evolving post-AGB stars 
can change the system’s UV color dramatically. 





Finally, most (but not all) stellar population codes neglect binary 
stars and binary evolution. In part, this is because modeling 
the flux distribution from a binary population involves lots of 


additional free parameters. But binary evolution probably is 


important! 
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Chemical Evolution 


Another topic related to the photometric evolution of a stellar 
population is the chemical evolution produced by a set of stars. 
The topic connects the evolution of single stars to the evolution 
of a galaxy. 





To start, let’s consider the types of parameters and variables that 
are involved. First, there are the global variables, all of which 
are a function of time. 





M,: ‘Total mass of interstellar gas 

M,: Total mass of stars 

Maw: Total mass of stellar remnants (white dwarfs) 
Mı: Total mass of the system 

E: the rate of mass ejection from stars 

Ez: the rate of metal ejection from stars 

W: the creation rate of stellar remnants. 


Naturally, M4 = Mg + Ms + Mw. 





Next, there are the global parameters of the model, which the 
investigator specifies. Again, all can be a function of time. 


W: Rate of star formation 

f: Rate of infall or outflow of material from the system 
Zr: Metal abundance of the infall (or outflow) material 
(m) : the Initial Mass Function 


As in the equations of photometric evolution, the IMF should be 
normalized so that the total mass is one, as in (31.02). 


Finally, you need four variables which come from stellar evolution 


w: the mass of a stellar remnant 

Tm: the main-sequence lifetime of a star 

Min: the turnoff mass of a population with t = T 

pz: the stellar recyclable mass fraction that is converted to 
metal z and then ejected into space. 





Given the above variables and parameters, the goal is to derive 
Z(t), the fraction of metals (individually, or as a group) in the 
interstellar medium as a function of time. 


Equations of Chemical Evolution 





There are five coupled differential equations which describe the 
chemical evolution of a system. 
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d(ZMy) = -ZÝ + Ez + Zf (32.05) 
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The equations are fairly simple to understand. Equation (32.01) 
is simple mass conservation. Equations (32.02), (32.03), and 
(32.04) keep track of the amount of mass that gets locked up 
in stars (or released into the ISM). Equation (32.05) is the most 
complex, as it describes how the metallicity of the interstellar 
medium changes with time. The first term of (32.05) refers to 
the amount of ISM metals that becomes locked up into stars, the 
second term gives the amount of metals being released by stars, 
and the third represents the amount of metals being brought in 
(or lost) from outside. 


Although the list of variables mentioned above is formidable, not 
all the variables are independent. Consider E, the mass ejection 
rate from stars. Since mass loss only occurs during post-main 
sequence evolution, the rate of mass loss is related to the number 
of stars evolving off the main sequence at any time. If the galaxy 





consisted only of a single population of stars, this rate would be 
given by 

dmin 
However, for a galaxy with on-going star formation, the calcula- 
tion of mass ejection rate must count the main-sequence turnoff 
stars of all stellar ages. Thus the total amount of ISM returned 
from stars at time t is 
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where M, is the upper mass limit of the stellar IMF, and Min, 
the turnoff mass at time t. Similarly, the equation for the total 
mass of remnants formed is 

Mu 


W = w(t — Tm) d(m,t — Tm)dm (32.07) 
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The equation for Ez is a bit more complicated since it has two 
terms: one to represent the amount of new metals created by 
a star and released during mass loss, and a second to represent 
the amount of metals that were lost from the ISM when the star 
formed, but are now being re-released. Mathematically, this is 
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Finally, there is an equation of metal conservation. If Z, is the 
average metal content in stars, then the total amount of metals 


produced in a galaxy over a Hubble time is 
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Primary and Secondary Elements 


The above equations assume that p,, the fraction of a star which 
is converted into metals, is independent of the initial metallicity 
of the star. In other words, it assumes that the metal under 
consideration is a primary element. However, some elements can 
only be made if another element already exists. For example, 
to make nitrogen (via the CNO cycle), the star must already 
have some carbon. Thus, the mass ejection rate of a secondary 
element, X, is 


Ex = i mpx Z(t — Tm)V(t — Tr) o(m, t — Tm)dm+ 
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Note that this is similar to the equation for Ez, in that it has 
two terms: the creation term and the recycle term. However, in 
this case, the creation term depends on the prior abundance of 
Z: 





Analytic Approximation to Chemical Evolution 


Obviously, solving the above coupled differential equations with 
their four free parameters (Y, ¢, f, and Zp) is a non-trivial nu- 
merical problem. However, the problem can be greatly simplified 
if you make two approximations. 


The first approximation to make is to say that the initial mass 
function of stars is independent of time. That is, (m,t) = o(m). 
Since little is known about how the IMF changes as a function of 
galactic conditions, this may, or may not, be a good assumption. 
(The prevailing theory is that the IMF for a system with no 
metals is heavily biased towards extremely massive stars, but 
once the first bits of metals get introduced into the ISM, this 
bias goes away.) 


The second approximation is call the instantaneous recycling ap- 
proximation and it is a bit tricker. The approximation says that 
there are two types of stars in a galaxy: those that live forever, 
and those that evolve and die instantaneously. Although this 
sounds like a poor assumption, it’s not as bad as it first appears. 
Recall that the timescales for stellar evolution: 





Main Sequence Lifetimes 


Spectral Mass Luminosity Lifetime 
Type (M/Mo) (L£/Lo) (years) 

O5 V 60 7.9 x 10° 5.5 x 10° 
BO V 18 5.2 x 104 9.4 x 10° 
B5 V 6 820 5.2 x 107 
AO V 3 54 3.9 x 108 
F0 V 1.5 6.5 1.8 x 10° 
G0 V 1.1 1.5 5.1 x 10° 
K0 V 0.8 0.42 1.4 x 101° 
MO V 0.5 0.077 4.8 x 101° 
M5 V 0.2 0.011 1.4 x 101 


Note the values. Stars with M > 5Mo evolve in less than 10° 
years, which, in cosmological terms, is almost instantaneously. 
On the other hand, stars with mass less than about 1Mo@ live 
forever. So the approximation only breaks down for a limited 
mass range. 


Let’s choose mı to be the dividing line between stars that live 
forever, and stars that evolve instantaneously. Let’s also define 
three new quantities, the Return fraction of gas 


R= f o COTS (32.11) 


the Baryonic Dark Matter fraction 


D= i wo(m)dm (32.12) 
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and the Net Yield (of element i) 
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In words, R is the amount of mass a generation of stars puts back 
into the ISM, D is the amount of mass a generation of stars turns 
into stellar remnants, and y; is the fraction of metal i produced by 
stars for every 1Mo of material locked up into stars or remnants. 
The importance of these three quantities is that each depend only 
on the IMF. If we assume some universal form for the IMF, then 
R, D, and y; are constants that depend only on stellar evolution. 
In other words, they are known quantities. 


Now, let’s take another look at the equations for E, W, and Ez. 
If we assume ¢(m) is independent of t and use the instantaneous 
recycling approximation, then 
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Similarly, the equation for stellar remnants becomes 
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~ v(t) if E w (32.15) 


and, after a bit of math, 


Ez =V{ZR+y,(1—R)} (32.16) 


With our two assumptions, the equations of chemical evolution 
become 
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a -(1-R-D)Y (32.18) 

Me = —(1—- R)Y¥ + f (32.19) 
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Actually (32.21) can be further simplified by noting that 


d(ZM4) „dM dZ 
dr = Zt Mg (32.22) 





Substituting (32.19) for dM, /dt then yields 


M, = = y,U(1— R)+ (Z; —Z)f (32.23) 


For secondary elements, the ejection rate from (32.10) becomes 


Ex =WVZ(1—R)y, + RXV —y.(1- R)XY 
= U(1— R)y,(Z —X)+UXR (32.24) 


(In many cases, X < Z, so to simplify things, you can often 
get away with Z — X ~ Z.) Metal conservation for secondary 
elements is therefore 

d(X M,) 
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But since 
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some of the terms cancel, so 


My = WI - R)(Z - Xue + (X - X)f 


(32.26) 


(32.27) 


Estimating the Stellar Yields 


Models of stellar evolution and supernova nucleosynthesis can 
predict the amount of metals produced by a stellar population. 
But is there a way to empirically check the results of these (purely 
theoretical) calculations? 


Let’s start with the equation for the total amount of metals pro- 
duced over a systems lifetime. Using the instantaneous recycling 
approximation, 


t Mau 
Z:(Ms + Mu) + ZM; = i J mp,U(t' — Tm)lm, t — Tm)dt'dm 
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= (1 — R) y-Yr (32.28) 


where V7 is the total amount of star formation over the history 
of the galaxy. Now that we’ve defined Ypy, let’s get rid of it. If we 
integrate the equations for the amount of mass that gets locked 
up in metals and remnants, 
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Now let’s define u as the gas fraction of a galaxy 
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So, if you choose a region (say, the solar neighborhood), and 
measure the average stellar metallicity, the present metal abun- 
dance in the interstellar medium, and the gas fraction, you can 
estimate the stellar yield of the elements in question. This is a 
key constraint on our understanding of the creation of metals by 
supernovae. 





Probing Star Formation 


The above equations allow us to estimate the history of star 
formation in the solar neighborhood. For example, if we re-write 
(32.18) as 

1 dMs 
1-R-D dt 


and split the derivative into three parts, we get 


1 dlog M, dlog Z dMs 
Į = — | ——— —— 2. 
1- R-D ( dlog Z ) ( dt ) (qe) eo) 


Now consider the derivatives. Since we can measure the metal- 
licity of solar neighborhood stars, we can determining how much 
stellar mass there is as a function of metallicity. Thus the first 
derivative in (32.35) is a measureable quantity. Similarly, if we 
study nearby F-stars, and compare their absolute luminosities to 
the F-star zero-age main sequence luminosity, we can estimate 
how much main-sequence evolution has occurred, t.e., we can es- 
timate their ages. If we measure the stars’ metallicities as well, 
then we have dlog Z/dt. Finally, the last derivative is simply 
In(10)M,. Thus, we can measure the history of star formation 
in the solar neighborhood. 


wv 





(32.34) 





This leads to the long-standing G-dwarf problem. Simple models 
of chemical evolution (say, a closed-box model or a constant in- 
fall model) predict many more low-metallicity stars in the solar 
neighborhood than are observed. 


To get the history of matter infall, we can do a similar manipu- 
lation with equation (32.19) 
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(32.36) 


If we measure the metallicities of clouds of H I and Hə of different 
masses, then all terms of this equation are known, and the history 
of matter infall can be found. 


The Closed Box Model of Chemical Evolution 


As an example of what a chemical evolution model can do, con- 
sider a closed system, where all the material for current star for- 
mation comes from mass lost by a previous generation of stars. 
In this case, there is no infall, and, from (32.23), 


dZ 


Mg =y.V(1—R)+(Z;-Z)f =y.¥(1-R) (32.37) 


In addition, from (32.19), we have 


a t= R)Y + f = -(1 — R)Y (32.38) 





By dividing these two equations, we get 
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Since y, is a constant of stellar evolution 
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(32.40) 
where Zp and M,, represent the initial metallicity and gas mass 
of the galaxy, and Z; and M,, represent those quantities today. 
Now let 
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or 
Zı — Zo 

Hı = Ho exp [-2-2) (32.43) 

In other words, as the system evolves, the gas fraction will de- 

crease exponentially with Z. Of course, we can’t measure the 

metallicity evolution of the ISM directly, but we can use stars as 

a probe. If we take the derivative of (32.43) with respect to Z, 
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Meanwhile, through (32.18) and (32.19) 
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Finally, if we put this equation in terms of log Z, instead of Z, 
then 


le mo (2) (852) (Bt) of 2-2 
dlog Z Ye 1—R O71 Yz 
(32.47) 
Because the number of stellar metallicity measurements has (tra- 
ditionally) not been extremely large, many times people plot the 
cumulative distribution, i.e., the number of stars with metallici- 


ties less than Z. This is simply found by integrating (32.47). If 
we collect the constant terms and let 


G= (=E) (=) (32.48) 








then for the cumulative distribution 
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When this is fit against observations, it is clear that the metal- 
licities of stars in the solar neighborhood cannot be fit with a 
closed box model of galactic evolution. Either f Æ 0, or Zo Æ 0, 
or the initial population of stars did not have the same IMF as 
the stars today, or there are severe chemical inhomogeneities in 
the ISM, and star formation occurs preferentially in regions with 
high metallicity. 





A similar closed-box calculation can be performed for secondary 
elements. For these, if you divide (32.27) by (32.19), then 
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If we assume that X < Z, then this simply reduces to 
Md X = —y,ZdM, (32.51) 


Now if we assume that M, = M; at t = 0, then from (32.42) 


Mga = M exp {=} (32.52) 


or 
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Thus 
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M,dX = ma = X=5 (| z? (32.54) 


In other words, if X is a secondary element, then a plot of log- 
X versus log-Z (i.e., |X] vs. [Z]) should have a slope of 2 (and 
presumably go through the solar value). 
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Fig 4.16 'Galaxies in the Universe’ Sparke/Gallagher CUP 2007 


Timescales for Chemical Evolution 


The timescales for chemical evolution are simple to derive. First, 
let’s consider the gas consumption timescale. In the absense of 
accretion, how long does it take to use up the gas? 


My My 
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In the solar neighborhood, Mg ~ 5.7Mo pe’, and the star 
formation rate is UV ~ 4.2Mo pe? Gyr’. For R ~ 0.4, that 
means that the gas should e-fold in 7, ~ 2 Gyr, much less than 
a Hubble time. 





Similarly, we can calculate the timescale for chemical enrichment. 
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Since yx ~ Z, Tz ~ Tą. In other words, for the solar neigh- 
borhood, the metallicity (and amount) of the gas should e-fold 
rather quickly. This strongly suggests that some our region of 
space has received matter from some other region. 


The Balanced Infall Model 


A closed-box model is not realistic for most systems: there is 
good evidence that infall of new material plays a part in the 
chemical evolution of a region. So let’s calculate (for simplicity) 
a model where the infall just balances the amount of material 
becoming locked up into stars. In other words, a model where 
the mass in the interstellar medium stays constant. In this case 


dt 





=-V+E+f=—-(1-R)U+f=0 (32.57) 


So 
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The equation for metals is then 


My = W(1— Ry. + (Zs — Z)f = fue + Zp- Z) (82.59) 


Also, since there is accretion 
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If we integrate this 
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Now let v represent the total amount of mass accreted, scaled to 
the mass in the ISM, i.e., v = (M —Mo)/M,. Then 


Z1 = (yz + Zi) — e”) + Zoe” (32.64) 
If the galaxy began with Zo ~ 0 and Zr ~ 0, then 
Zi =yz(1-—e ”) (32.65) 


Note that if you also assume that the galaxy began as an entirely 
gaseous system, then 
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Thus, as u > 0, v > œ, and Zı — yz. In other words, the 
metallicity of the system asymptotes out at the value of the stellar 
yield. 


The balanced infall model is similar to the closed-box model in 
that Z x yz, so the metallicity of a system is actually measuring 
the stellar yields. Also, note that Z/y, is not a strong function 
of the gas fraction, or Y/f. It depends mostly on the current 
properties of the system. 


